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The  commonly  accepted  linear-parabolic  oxidation  model  for  the  thermal  oxidation  of 
Si  includes  two  rate  processes  in  a  steady  state:  reaction  between  Si  and  oxidant  at  the 
Si-SiO^  interface;  and  transport  of  oxidant  through  the  SiOj  film.  Based  on  available 
data,  it  is  argued  that  the  former  process  seems  dominant  for  thin  film  growth  in  dry  0^. 
A  number  of  measured  SiC^  film  and  Si-SiC^  interface  measured  properties  are  reported,  as 
well  as  the  variation  of  the  these  properties  with  oxidation  temperature  and  Si  substrate 
orientation.  These  properties  include:  refractive  index,  density,  intrinsic  stress, 
interface  fixed  oxide  charge  and  interface  trapped  charge.  It  is  also  observed  that  all 
of  these  properties  display  similar  oxidation  temperature  and  inert  anneal  behavior  plus 
a  complex  orientation  dependence.  Through  the  use  of  a  modified  form  for  the  interface 
reaction,  a  better  understanding  is  obtained  of  both  the  origin  of  these  measured 
properties,  and  new  oxidation  data  taken  on  five  Si  orientations  and  at  lower  oxidation 
temperatures  than  previously  reported. 


The  study  of  the  thermal  oxidation  of  single  crystal  silicon  remains  an  important 
area  of  research  within  the  field  of  electronic  materials.  This  is  so  because  of  both  the 
technological  relevance  of  the  Si-SiO^  interface  (1,2)  with  the  associated  important 
scientific  questions  pertaining  to  oxide  film  growth  kinetics  (3,4),  and  the  general 
questions  pertaining  to  the  origin  and  control  of  electronics  properties  of  semiconductor 
surfaces . 

The  present  paper  first  compares  the  relative  role  of  the  interface  reaction  between 

Si  and  oxidant  at  the  Si-SiO^  interface  with  the  transport  of  oxidant  by  steady  state 

Fickian  diffusion  to  the  interface.  Using  available  diffusivity  values,  D,  parabolic  rate 

constants,  k  and  oxidation  rate  data,  L,  (5-7),  it  is  found  that  the  interface  reaction 
P 

'for  Si  oxidation  is  dominant  up  to  at  least  several  tens  of  nm  of  SiO^  film  growth  with 
oxidation  temperatures  from  600°C  up  to  1000°C  in  1  atm  dry  0 Secondly,  we  report  the 
measurement  of  a  number  of  physical  properties  for  SiO^  films  on  Si  substrates:  SiO^  film 
density,  refractive  index,  intrinsic  film  stress  and  fixed  oxide  and  interface  trapped 
charge.  Some  of  these  measurements  are  recently  reported,  others  have  been  in  the 
literature.  The  oxidation  temperature  dependence  for  all  of  these  properties  is 
qualitatively  the  same,  and  strongly  suggests  a  common  origin.  Thirdly,  we  consider  the 
interface  reaction  between  Si  and  oxidant  and  deduce  a  model  for  this  reaction  to  explain 
initial  oxidation  kinetics.  New  oxidation  data  (6,7)  is  presented  that  shows  some 
agreement  with  the  model  which  explicitly  considers  both  the  areal  density  of  Si  atoms  on 
various  Si  surfaces  and  the  intrinsic  stress.  The  stress  is  compressive  in  SiO^,  thereby 
reducing  the  diffusivity  of  oxidant  through  the  SiO^  film;  and  the  stres  is  tensile  in  Si, 
thereby  enhancing  the  reactivity  of  the  Si  surface,  if  there  is  any  effect  of  stress  at 
all.  Experimental  detail  is  treated  in  separate  publications;  only  a  brief  summary  of 
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relevant  experimental  results  are  presented  herein. 


I.  The  Initial  Regime 

There  are  numerous  publications  (8-10)  which  attempt  to  consider  the  transport  of 

oxidant  to  be  always  rate  limiting  for  the  oxidation  of  Si.  We  find,  however,  that  this 

assumption  is  not  justified  and  we  base  our  conclusion  on  a  comparison  of  the  calculated 

diffusion  flux  with  the  actual  oxidation  rate.  Using  available  data  for  D,  k  and  L  and 

P 

the  assumption  of  steady  state  Fickian  diffusion,  the  diffusion  flux  of  0^,  F(D),  is  given 
by: 


F(D)  =  Dq  (Cj-C^/L 

where  is  the  solubility  of  0^  in  SiO^  and  is  the  0^  concentration  at  the  Si-SiO^ 

interface.  We  assume  that  the  solubility  of  0^  in  SiO^  is  constant  with  temperature  and 

16  —  3 

has  the  value  Cj  =  5  X  10  cm  (11).  We  need  to  realize  that  the  D  values  of  Norton  (5) 
were  obtained  at  high  temperatures  and  using  bulk  SiO^  glass,  and  thus  cannot  be  readily 
extrapolated  to  thin  films  of  pure  SiC^  prepared  at  temperatures  below  900°C  where  the 
intrinsic  compressive  6tress  (12,13)  may  seriously  decrease  D  (8-10).  For  these  lower 
temperatures,  we  calculate  D  from  the  available  k^  data  (6)  and  with  the  use  of  the  Deal 
and  Grove(ll)  expression  for  D  as: 


D  =  Nk  /2C, 

P  1 

where  N  is  the  number  density  of  *n  SiO^.  This  method  of  calculating  D  values  for 

lower  temperature  oxidation  is  in  accord  with  the  experimentally  obtained  curvature  of  the 
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In  vs.  1/T  plot(14)  which  shows  a  negative  deviation  to  the  linearity  extrapolated  from 
high  temperatures  to  the  lower  temperatures  (below  900°C).  This  results  in  a  smaller  value 
for  D  than  the  values  extrapolated  from  Nortons  data.  We  assume  is  very  small  relative 
to  Cj  and  can  therefore  be  ignored. 

We  compare  this  diffusive  flux,  F(D),  with  the  flux  of  0^  calculated  using  the 
experimental  Si02  growth  rates,  F(exp).  A  ratio,  R,  is  formed: 

R  =  F(D)/F(exp) 

This  ratio  has  the  following  characteristics: 

R  =  1  Diffusion  control 
R  >  1  Interface  control 
R  <  1  Impossible 

and  can  therefore  be  used  to  indicate  the  dominant  kinetic  process.  When  the  diffusion 
flux  fully  accounts  for  the  actual  film  growth  rate  then  diffusion  is  the  rate  limiting 
process  and  R  =  1.  When  the  calculated  diffusion  flux  is  larger  than  is  necessary  to 
account  for  the  measured  oxidation  rate,  L,  then  the  oxidation  mechanism  is  under 
interface  reaction  control  and  R  >  1.  The  diffusive  supply  function  can  never  yield  less 
C>2  than  is  required  by  L  hence  R  <  1  is  not  possible  for  our  situation  as  defined.  The 
results  of  this  calculation  are  displayed  in  Table  1.  It  is  seen  that  R  is  always  greater 
than  unity  for  oxidation  temperatures  of  from  600°C  to  1000°C  and  up  to  at  least  several 
tens  of  nm  which  comprises  the  initial  fast  oxidation  regime.  Therefore,  based  on  the 
available  data  and  simple  ideas  of  transport  and  reaction,  we  conclude  that  the  reaction 


at  the  Si  -  SiC^  interface  is  likely  to  be  rate  limiting.  There  exists  other  supporting 
evidence  for  this  position.  For  example,  the  oxide  intrinsic  stress  is  compressive  in 
SiOj  and  tensile  in  Si.  Thus  if  there  is  any  effect  of  stress  for  thin  films,  the  stress 
should  decrease  the  diffusivity  by  decreasing  the  free  volume  in  the  Si02(8-10),  and 
increase  the  surface  reactivity  by  stretching  the  Si  bonds.  We  actually  observe  an 
enhanced  initial  oxidation  rate  relative  to  the  kinetics  for  thick  films  which  seems  to 
support  the  interface  reaction  effects.  Also,  if  the  diffusivity  is  rate  limiting  and 
affected  by  stress,  then  in  the  early  stage  of  oxidation  the  diffusivity  must  be  smaller 
than  for  thick  films  where  some  oxide  has  vi scoel as t ica 1 1 y  relaxed,  yet  a  larger  initial 
rate  is  observed  for  dry  0  oxidation. 


II.  Physical  Properties  Measurements 

We  have  measured  several  physical  properties  of  SiC>2  films  on  Si  and  there  exists  a 
number  of  similar  measurements  in  the  literature.  We  select  for  display  in  this  paper 
only  those  properties  whose  values  we  have  confirmed,  and  more  importantly,  we  show  the 
variation  of  these  properties  with  oxidation  temperature. 

Figure  1  shows  the  measurement  of  the  real  part  of  the  refractive  index  for  Si02  as 
function  of  oxidation  temperature.  This  was  assumed  by  Taft  to  be  a  result  of  an 
increased  oxide  density  (15).  Later  measurement  of  the  index  and  density  (16)  revealed 
almost  exactly  the  same  values  for  the  refractive  index,  n,  as  Taft.  While  the  density 
measurement  was  not  nearly  as  accurate,  a  definite  increase  was  found  in  the  measured 


density  for  the  lower  oxidation  temperatures.  Figure  2  shows  the  results  of  converting 

p,  through  the  Lorentz  -  Lorentz  equation: 


p  =  K( n2- 1 ) / ( n2+2 ) 


the  n  values  to  densities, 


and  in  order  to  calculate  K,  literature  values  for  n  at  measured  p  were  used  (17)  from 
which  a  value  of  K  *  8.0461  was  obtained. 

Further  evidence  for  densif icat ion  was  recently  obtained  using  infrared  spectroscopy 
techniques  (18).  Figure  3  shows  a  typical  set  of  spectra  from  which  the  dominant  change 

is  a  shift  in  the  1077cm  *  band  to  lower  frequencies  with  decreasing  oxidation 

temperatures.  This  band  is  attributed  to  the  in  plane  stretching  of  Si-O-Si  and  the  shift 
has  been  attributed  to  a  decrease  in  the  Si-0-Si  bond  angle  which  is  anticipated  for 
densif icat ion  (19). 

Figure  4  shows  the  results  of  intrinsic  film  stress  measurements  on  SiO^  films  on  Si 
(20,21).  The  details  of  the  actual  stress  measurement  technique  and  thermal  stress 
corrections  were  published  separately  along  with  the  experimental  results  (20,21).  These 
results  agree  essentially  with  previous  independent  measurements  (12,13).  We  observe  an 
increase  in  the  compressive  SiO^  intrinsic  stress  with  decreasing  oxidation  temperatures 
for  all  orientations.  For  reasons  to  be  discussed  later,  the  (111)  yields  the  smallest 
stress,  although  this  orientation  is  anticipated  to  yield  the  largest  stress  on  a  planar 

Si  surface.  The  film  stress,  o~ ,  in  the  elastic  limit  is  given  as: 

O”  —  E4 

where  E  is  Young's  modulus  and  4  is  the  strain  which  is  the  change  in  volume  divided  by 
the  volume,  oV/V.  The  orientation  dependence  of  o-  is  obtained  from  only  the  orientation 
dependence  of  E  (or  better  E/(l-v)  where  v  is  Poisson's  ratio),  since  the  volume  change, 
oV,  due  to  oxidation  is  the  same  for  all  orientations.  Table  2  shows  the  calculation  of 
E/(l-v)  values  for  various  Si  orientations  (7).  From  this  table  and  the  above 


relationship,  we  deduce  that  the  order  for  the  stress  on  the  various  Si  orientations 
should  be: 

°~1 1 1  >  °"l  10  >  “"lOO* 

However,  we  measure  the  order  for  the  stresses  to  be: 

°~1 10  ^  100  >  ^Ill- 

Later  we  discuss  both  the  possible  reasons  for  this  order  change,  and  more  importantly  for 
our  present  purpose  are  the  implications  of  this  experimentally  determined  orientation 
dependence  of  the  intrinsic  stress  on  Si  oxidation  kinetics. 

The  remaining  physical  properties  reported  here  include  the  fixed  oxide  charge,  Qp 
and  interface  trapped  charge,  Q^.  The  former  charge,  Q^,  is  associated  with  excess  Si  at 
the  Si-SiO^  interface  while  is  the  charge  trapped  in  surface  electronic  states.  These 

surface  states  are  thought  to  be  a  result  of  unsatisfied  or  dangling  bonds  at  the  Si 
surface.  Fig.  5  is  a  pictorial  representation  of  the  trend  observed  for  and  Q.^  with 
oxidation  temperature  (22,23).  Both  of  these  kinds  of  charges  have  been  found  to  anneal 
to  lower  values,  but  quite  complex  chemical  effects  are  sometimes  observed  for  Q.  . 

It  is  important  to  realize  that  all  the  properties  thus  far  discussed  decrease  in 
value  with  inert  gas  annealing  (16,21-23).  Also,  all  the  property  values  are  lower  for 
oxides  grown  in  an  H^O  containing  ambient.  The  rate  of  decrease  of  the  property  seems  to 
depend  on  the  temperature  of  the  anneal.  It  therefore  appears  that  there  is  a  relaxation 
upon  annealing  which  is  strongly  suggestive  of  viscoelastic  behavior ( 12 , 13 , 13 ) . 


Fig.  6  summarizes  all  the  herein  reported  physical  property  results.  With  this 
similar  behavior  of  the  properties  with  oxidation  behavior,  we  believe  that  a  common 
origin  is  strongly  suggested.  We  consider  that  the  stress  is  the  most  fundamentally 
related  property  to  the  oxidation  reaction  through  the  change  in  the  molar  volume,  ^V. 

This  change  in  volume  causes  the  SiC^  to  increase  in  density  through  two  mechanisms.  One 
is  a  simple  compression  due  to  the  existence  of  a  residual  stress  as  measured.  This 
effect  is  quantified  from  the  stress  optical  coefficient  for  SiO^  (24). 

dn/do-  *  9  x  10  ^  cm^/dyne 

From  this  coefficient,  we  can  only  account  for  less  than  half  of  the  observed  change  in  n. 
The  remaining  dens i f icat ion  may  be  due  to  structural  changes,  perhaps  an  alteration  of  the 
ring  statistics  (25).  This  latter  effect  may  account  for  the  lack  of  complete  annealing 
of  some  of  the  propert ies ( 21  ) .  The  rationale  used  here  is  that  the  stress  induced 
components  that  require  the  continual  existence  of  the  stress  are  annealed  as  the  stress 
is  annealed,  but  the  structural  changes  that  required  stress  initially  and  for  relaxation 
require  bond  breaking,  are  less  likely  to  anneal  quickly  or  completely. 

With  Part  I  indicating  the  importance  of  the  interface  reaction  and  Part  II 
demonstrating  that  a  number  of  properties  related  to  the  film  stress  change  with 
decreasing  oxidation  temperature,  we  now  use  the  effect  of  stress  on  the  interface 
reaction  between  oxidant  and  Si  to  attempt  to  correlate  previously  published  and  newly 
obtained  oxidation  data. 

III.  Oxidation  Kinetics 
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In  this  part,  we  firstly  present  oxidation  data  of  Massoud  et  al .  (6)  which  shows  a 


crossover  effect  with  various  Si  orientations.  Fig.  7a  and  b  shows  that  initially  the 
(110)  orientation  has  the  greatest  rate  of  oxidation,  but  at  thickness  greater  than  20  run 
the  (111)  is  dominant.  Previously,  we  attempted  to  explain  this  phenomenon  (25)  based  on 
a  recent  revision  (26)  to  the  linear  -  parabolic  model  (11).  Only  the  essential  features 
of  the  revised  model  as  it  pertains  to  the  initial  regime  will  be  reproduced  here.  Fig. 
shows  that  as  the  oxide  grows,  the  molar  volume  change  and  the  adhesion  to  the  Si 
surface  cause  the  development  of  a  compressive  oxide  stress  (tensile  in  the  Si  surface) 
and  a  viscous  flow  in  the  normal  direction  which  is  parallel  to  the  SiO^  growth  direction 
Using  this  picture,  and  starting  from  the  Deal-Grove  (11)  interface  reaction  flux: 


¥ 


klC2 


where  is 
the  Si-SiO^ 
of  Si  atoms 


first  order  rection  rate  constant  and  is  the  concentration  of  oxidant  at 

interface,  we  modify  this  flux  by  explicitly  expressing  the  available  number 

★ 

on  the  Si  surface,  C  .  ,  as: 

O  1. 


F 


C2CSi 


★ 


* 

atlC*  ^Si  *s  °^ta^ne^  directly  from  the  geometric  areal  density  of  Si  atoms  on  a  given 
surface  ,  CSl’  by  multiplying  Cg^  by  a  fraction.  Considering  SiC>2  as  a  Maxwell  solid,  we 
assume  that  this  fraction  is  the  strain  rate,  4  as: 


where  cr  is  the  stress  and  the  oxide  viscosity.  This  means  that  the  rate  at  which  the 
SiC^  relaxes  into  the  growth  direction  influences  the  rate  at  which  the  oxidation  reaction 
will  proceed.  Then  the  following  form  for  from  the  Deal-Grove  expression  is  obtained: 

ki  =  k'csiwr! 

It  should  be  realized  that  this  modified  expression  for  k^  is  obtained  by  simply  assuming 
that  stress  affects  the  interface  reaction  and  that  the  number  of  Si  atoms  on  the  surface 
is  also  infuencing  the  oxidation  rate  in  the  early  regime.  Even  if  our  assumptions  are 
entirely  correct,  it  is  quite  likely  that  the  form  for  the  rate  expression  will  require 
modification  as  more  data  is  revealed.  However,  this  simple  formulation  will  enable  a 
scaling  with  the  oxidation  data  and  in  particular  the  crossover  effect  shown  above. 

In  order  to  explain  the  Massoud  et  al .  (6)  data  using  the  relationship  above,  we 
assume  that  in  the  earliest  stage  of  oxidation,  the  primary  factor  is  the  areal  density  of 
Si  atoms.  As  seen  in  Table  3,  the  oxidation  order  from  this  is: 

(110)  >  (111)  >  (100) 

and  this  scales  with  the  data.  As  the  oxide  film  grows,  the  third  dimension  of  the  SiO^ 
network  develops,  and  with  this  we  should  anticipate  the  film  to  exert  itself 
mechanically.  Using  the  revised  form  for  kj,  we  would  then  expect  the  stress  to  dominate. 
At  the  time  we  submitted  this  work  for  pub  1 ica t ion ( 26 ) ,  we  did  not  haveaccess  to  the 
measured  stress  values  in  Figure  4,  and  thus  we  relied  on  the  calculated  stress,  thereby 
anticipating  that  the  order  of  the  oxidation  kinetics  based  on  stress  should  be: 


(Ill)  >  (110)  >  (100) 


which  is  as  observed  for  the  oxide  grown  thicker  than  the  crossover.  Of  course,  now  with 
the  stress  measurements  as  a  function  of  Si  orientation  (Figure  4),  we  find  that  this 
order  is  not  correct.  Also,  we  have  now  added  two  new  orientations,  the  (311)  and  (511), 
which  according  to  Table  3  and  the  atom  counting  procedure  of  Ligenza  (28)  should  have  a 
greater  areal  density  of  Si  and  hence  a  greater  initial  oxidation  rate.  This  new 
oxidation  rate  data  which  includes  these  additional  Si  orientations  is  shown  in  Figure  9. 
We  observe  that  the  (511)  and  (311)  rates  are  respectively  further  above  the  (100)  rate, 
but  both  are  below  the  (111)  rate.  Ligenza  presumably  considered  the  (511)  and  (311)  to 
be  simple  atomic  planes  and  counted  the  atoms  on  or  near  the  plane.  In  fact,  however,  the 
(511)  and  (311)  planes  are  vicinal  planes  to  the  (100)  and  (111)  (29).  As  such,  both  of 
these  planes  are  composed  of  a  combination  of  (100)  terraces  and  (111)  ledges  in  ratios  so 
as  to  yield  the  proper  inclination  to  the  (100)  plane.  The  proportion  of  (111)  steps 
increases  from  (511)  to  (311)  to  eventually  produce  the  (111)  and  the  oxidation  rate 
scales  similarly.  Thus  we  continue  to  believe  that  the  initial  oxidation  rate  scales  with 
the  areal  density  of  atoms,  and  that  apparently  an  error  exists  in  the  literature  for  the 
areal  density  for  the  (311)  plane. 

According  to  the  revised  interface  reaction  expression,  the  stress  at  the  Si  surface 
should  also  be  important.  Indeed,  the  product  Cg^o-  ought  to  be  determined  through  the 
linear  rate  constant.  However,  Table  4  shows  that  this  product  is  anomolously  low  for  the 
(111)  surface.  Since  the  very  initial  rate  scales  with  Cg^,  we  continue  to  support  the 
previous  assumption  that  the  stress  in  the  very  initial  stage  is  not  important.  We 
previously  argued  (26)  that  this  was  so  because  the  film  was  too  thin  to  mechanically 
exert  itself  until  the  third  dimension  grows  adequately  thick  and  hence  to  develop 


.  •  .  ' .  •  ,'r  ■-  •-  .v  .v  - 


mechanical  properties.  Other  possible  explanations  are  that  the  stress  for  very  thin 

films  is  different  than  for  the  films  of  greater  than  100  nm  as  in  Figure  4.  We  have 

commenced  experiments  directed  at  films  less  than  20  nm.  Also,  as  stated  above,  the  form 

for  the  revised  rate  constant  may  not  be  correct.  Another  possibility  is  that  the  low 

stress  for  the  (111)  surface  may  be  responsible  for  an  enhanced  rate  due  to  an  enhanced 

diffusivity  rather  than  an  effect  on  the  surface  reaction  kinetics.  Evidence  for  this  is 

given  in  Table  5  which  contains  a  set  of  parabolic  rate  constants  from  Massoud  et  al.  (6). 

These  constants  show  that  for  lower  temperatures,  the  k  values  for  the  (111)  actually 

P 

increase  above  those  for  the  (110)  and  (100)  orientations  thus  indicating  the  importance 
of  transport  but  somehow  orientation  dependent  transport.  Of  course,  since  these  k 

P 

values  originate  with  the  Deal-Grove  model,  they  are  model  dependent  parameters  and  are 
therefore  only  as  reliable  as  the  model. 

Lastly,  we  return  to  the  stress  measurements  and  the  possible  reasons  why  the  (111) 
exhibits  a  lower  stress  than  anticipated  from  E/l-v  values  as  shown  in  Table  3.  The 
anticipated  stress  values  assumed  a  planar  Si.  surface.  Mott  (30)  has  suggested  that  if 
oxidation  proceeds  at  surface  steps,  the  resulting  stress  would  be  considerably  reduced 
because  a  part  of  the  is  relieved  by  the  receeding  Si  step  as  shown  in  Figure  10.  Hahn 
and  Henzler  (31)  report  that  among  a  rather  wide  list  of  process  parameters  that  alter  the 
surface  step  density  on  Si  surfaces,  they  measure  a  significantly  greater  density  of  steps 
as  the  (111)  compared  to  the  (100)  surface.  This  finding  lends  considerable  support  to  the 
application  of  the  model  proposed  by  Mott.  Contained  within  the  papers  at  this  workshop 
ts  one  of  particular  application  to  the  present  research  presented  by  Leroy  (32)  which 
quantitatively  explains  the  experimentally  found  orientation  dependence  of  the  stress. 
While  we  do  not  reproduce  the  Leroy  arguments  here  as  the  paper  is  published 
simultaneously  with  the  present  paper,  a  brief  summary  is  in  order.  The  Leroy  model 


considers  that  the  first  oxidation  event  creates  a  step  on  the  Si  surface,  since  Si  on  the 
surface  is  replaced  with  SiO^.  Once  a  step  is  produced,  the  forces  resulting  from  the 
oxidation  induced  will  be  manifest  in  the  residual  measured  intrinsic  stress  according 
to  the  resolved  component  of  the  biaxial  stress  in  the  Si  surface.  This  stress  requires 
resolution  in  the  direction  of  the  slip  plane  in  relation  to  its  angle  with  the  Si 
surface.  This  angle  is  orientation  dependent.  The  resulting  stress,  cr  is  compared  to  the 
critical  shear  stress,  X^.  When  o-  exceeds  X  ,  defects  are  produced  and  stress  is 
relieved.  Leroy  calculates,  quantitatively,  the  relative  order  for  the  (100),  (110), 

(111)  and  (311)  stresses  using  these  ratios  for  the  different  Si  orientations  of  this 
study.  He  considers  that  any  defects  produced  are  rapidly  removed  by  oxidation.  There 
are  further  profound  implications  of  the  Leroy  n»del  and  the  reader  is  referred  to  the 
paper  (32). 

Summary  and  Conclusions 

Using  available  data  for  D,  and  L  and  simple  calculations,  we  conclude  that  the 
interface  reaction  between  Si  and  oxidant  is  important  for  the  early  thermal  oxidation 
regime.  Other  evidence  such  as  possible  stress  effects  and  the  apparent  non-physical 
results  obtained  by  using  the  interface  reaction  further  strengthen  this  approach.  A 
variety  of  SiC^  and  Si-SiO^  physical  properties  have  very  similar  oxidation  temperature 
behavior  and  the  property  values  relax  to  lower  values  upon  inert  anneal.  This  suggests 
that  the  properties  have  a  common  origin  and  we  believe  that  the  development  of  the  stress 
is  the  common  cause  for  the  temperarture  dependence  of  the  properties.  The  anneal  behavior 
suggests  the  viscoelastic  nature  of  the  Si02  films.  A  revised  linear  rate  constant  is 
employed  along  with  new  oxidation  data  that  includes  five  Si  orientations  to  understand 
the  low  temperature  oxidation  phenomena.  The  importance  of  the  areal  density  of  Si  atoms 
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appears  to  be  a  strong  correlation  for  the  very  initial  regime  and  is  contained  explicitly 
in  the  revised  interface  model.  The  observed  crossover  seems  to  be  due  to  stress  transport 
effects  if  the  data  is  explained  within  the  proposed  revised  model.  While  the  revised 
interfacial  reaction  model  has  some  physical  appeal,  more  research  is  needed  to  prove  or 
disprove  this  model  and  to  arrive  at  a  better  description  of  the  early  stage  of  thermal 
oxidation  of  Si. 
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